Room temperature ballistic transport in narrow graphene strips 
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We investigate electron-phonon couplings, scattering rates, and mean free paths in zigzag-edge 
graphene strips with widths of the order of 10 nm. Our calculations for these graphene nanostrips 
(GNSs) show both the expected similarity with single-wall carbon nanotubes (SWNTs) and sup- 
pression of the electron-phonon scattering due to a Dirichlet boundary condition that prohibits one 
major backscattering channel present in SWNTs. Low-energy acoustic phonon scattering is expo- 
nentially small at room temperature due to the large phonon wavevector required for backscattering. 
We find within our model that the electron-phonon mean free path is proportional to the width of 
the nanostrip and is approximately 70 /im for an 11-nm-wide nanostrip. 

PACS numbers: 73.50.-h, 73.23.Ad, 73.61.Wp 
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I. INTRODUCTION 

Low-dimensional systems have long attracted attention 
in mesoscopic physics. One area of such interest is one- 
dimensional conductors. These quantum wires do not 
obey the usual macroscopic Ohm's Law. Instead they ex- 
hibit ballistic or coherent diffusive electron transport. In 
the former transport regime, the conductance is indepen- 
dent of the length of the quantum wire and is quantized 
in units of 2e 2 /h. The requirement for ballistic transport 
is that the length of the quantum wire is shorter than 
its characteristic mean free path. In metallic SWNTs 
at room temperature, low-energy mean free paths of or- 
der one micrometer have been observed 1 ^. These ob- 
servations are also confirmed by theory, which assumes 
that electron-phonon scattering is the dominant source of 
scattering^ as the effects of static disorder are strongly 
suppressed 3 -^! 6 -. 

In this paper, we present electron-phonon calculations 
of graphene nanostrips terminated by hydrogen atoms. 
The interior of the nanostrips consists of an sp 2 -bonded 
graphene honeycomb lattice, and therefore there is hope 
that the nanostrips will inherit some of the special prop- 
erties of graphene. There are several studies of the elec- 
tron properties of GNSsLS£^lL12^ii£^L12^2i2a, 
and recent experimental progress 21 - in fabricating the ma- 
terials will likely add further interest. The nanostrips can 
be made from a sheet of graphene using lithography. Un- 
til recently, it was uncertain whether graphene would be 
thermodynamically stable or would spontaneously curl 
up into scrolls 22,23 . Initial electron transport measure- 
ments on graphene have demonstrated high room tem- 
perature mobilities, suggesting that a sheet of graphene 
is relatively inert to interface scattering at the sub- 
strate surfac o 23 ! 24 ' 25 . From experiments, the graphene 
mean free path has been estimated to be 400 nm 2 ^ and 
600 nm 21 . The mean free path has not shown significant 
temperature dependence, and therefore we would expect 
the electron-phonon mean free path to be even longer. 

Unlike in graphene, the carriers in GNSs are confined 
to one dimension. The electron energy dispersion is quan- 



tized in the transverse direction like the dispersion in car- 
bon nanotubes. However, there are two important differ- 
ences between carbon nanotubes and nanostrips. First, 
metallic carbon nanotubes exhibit two channels (exclud- 
ing the spin degree of freedom) at the Fermi level due to a 
periodic boundary condition in the transverse direction. 
Nanostrips, on the other hand, have a Dirichlet boundary 
condition and only a single channel. Second, zigzag-edge 
nanostrips have edge states which do not exist in carbon 
nanotubes^. These edge states have energies close to 
the Fermi level, ep, and can therefore significantly affect 
low-bias electron transport. On the other hand, the edge 
states are highly localized on the edge o 7 ' 8 ' 10 and have 
only small overlaps with the extended states, thus leav- 
ing the extended states in the dispersion relatively intact. 
Similar edge localization has been seen in phonon energy 
dispersions 26 . An important consequence of the presence 
of edge states is that they replace low-energy channels 
which would otherwise cause backscattering. 

The next section describes the electron-phonon cross- 
sections in our model. These results are then used in 
Sec. lIIII where Fermi's Golden Rule is applied to calcu- 
late scattering rates and the electron-phonon mean free 
path. Our findings are finally discussed in the concluding 
section. 



II. ELECTRON-PHONON COUPLINGS 

We model the electron using the usual tight-binding 
Hamiltonian 
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where the summation is over all nearest neighbors 
■fij are hopping parameters, and c] and c; arc the 
fcrmionic creation and annihilation operators, respec- 
tively, on site i. When there is no lattice deforma- 
tion, we use the graphene 7r-orbital hopping parameter 
lij = 7o ~ — 2.7eV— . The energy dispersion for a nano- 
strip with N zz = 50 interior zigzag lines along the strip 
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FIG. 1: Electron energy dispersion of a graphene nanostrip 
(solid bands). This particular nanostrip has a width of N zz = 
50 interior zigzag lines which corresponds to approximately 
11 nm. The dispersion is compared to that of a (25, 25) carbon 
nanotube (dashed bands). Due to the edge states there is only 
one channel close to the Fermi level which is located at E = 0. 
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FIG. 2: First Brillouin zone with (a) quantized electronic 
bands and (b) phonon branches satisfying q x = 0. Phonons 
at the K-point can scatter an electron in state K into K'. 



axis is shown in Fig.[TJ Close to the Fermi level, the 
strip has only a single channel in each direction, and the 
channels are separated by a wavevector approximately 
q = (4/3)(7r/a), where a 0.246 nm is the lattice spac- 
ing. As a consequence, long- wavelength acoustic phonons 
cannot cause backscattering in zigzag-edge nanostrips. 

In this paper, we restrict ourselves to single-phonon 
events with first-order variations in the hopping pa- 
rameters. We limit the calculation to the energy gap 
AE w (25/N zz ) eV in which there is only a single channel 
in each direction. In graphene with a transverse bound- 
ary condition, there are only certain allowed transverse 
wavevectors as indicated in Fig.O Only phonons with 
the same set of wavevectors can efficiently scatter elec- 
trons due to conservation of crystal momentum. Fur- 
thermore, low-energy intraband electron backscattering 
can only be caused by phonons in the proximity of the 
K-point. We will only consider the case where q x = 
as all other equivalent K-points can be translated onto 
this band by a reciprocal lattice vector. Because of the 
edge states in zigzag-edge nanostrips, the quantization 
of the two-dimensional Brillouin zone is no longer exact 
although it still remains a good approximation in wide 



FIG. 3: Lattice vibrations which cause non-negligible intra- 
band electron-phonon coupling. The two modes represent (a) 
a LA phonon and (b) a TO phonon at a K-point. Note that 
due to the LA and LO degeneracy, the LA mode has been 
obtained in the limit q —> (4/3) (w /a). 



nanostrips for |fc| < (2/3)(7r/a), as the edge states are pri- 
marily located outside this wavevector windowZ^. There 
are six phonon branches with q x = 0. In two of these 
branches, the phonon modes arc entirely out-of-plane. 
These modes do not affect the bond length to first or- 
der and can therefore be neglected. The remaining four 
phonon modes are longitudinal acoustic (LA), longitu- 
dinal optical (LO), transverse acoustic (TA), and trans- 
verse optical (TO) phonons. The LA and TO phonons 
cause intraband scattering and the LO and TA phonons 
cause intcrband scattering between the conduction and 
valence bands. As we will see below, the interband scat- 
tering can be minimized by doping. This leaves the LA 
and TO branches with corresponding modes at a K-point 
shown in Fig.[3l 

We can estimate the electron-phonon coupling in the 
nanostrips by introducing a frozen phonon into the lat- 
tice. The frozen phonon perturbs the lattice and causes 
local variations in the bond lengths, which in turn affect 
the nearest-neighbor hopping parameters. We model this 
effect by expanding the Goodwin-Skinner-Pettifor scaling 
function for carbor>2£ and obtain hopping parameters of 
the form 



Jij = 7o - a dij ■ (ui - Uj), 



(2) 



where the coupling constant is a w 52eV/nm, is 
the unit vector in the direction to the nearest neigh- 
bor, and Ui — Uj is the relative lattice displacement. An 
LA or TO phonon at the K-point opens up a gap at 
k = ±(2/3)(7r/a). This effect can be seen in the reduced 
zone dispersion in Fig.[DJa), where k = ±(2/3)(7r/a) has 
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on an extended state. The figure also demonstrates that 
the conduction band is similar to that in a nanotube for 
|fc| < (2/3)(n/a) which is also expected as the eigenstates 
become identical in the limit N zz — ► oo. The transition 
between the edge and extended states at \k\ (2/3)(7r/a) 
also become sharper for wider nanostrips. Motivated by 
Fig-HJc) we take the electron-phonon matrix elements in 
graphene as an upper bound of those in GNSs and will 
from here on use the former, which we can estimate an- 
alytically. 

The electron-phonon interaction in Fourier space can 
be expressed as 



cl— ph 



[ikq 



(3) 



where M% q are the coupling matrix elements and 4>^ is 
the phonon operator of mode [i. The latter is given by 

1/2 

(4) 



where a'i, a_ q are bosonic creation and annihilation op- 
erators, and ujff are phonon frequencies. The phonon 
dispersion is calculated numerically using a force con- 
stant model^ with new parameters^. Following the 
derivations of matrix elements in armchair carbon 
nanotube o 30 ' 31 we find the intraband matrix elements ap- 
propriate for nanostrips 



FIG. 4: Effects due to a frozen TO phonon in a GNS 
[Nzz = 50) lattice, (a) The electron dispersion of the nano- 
strip is calculated with the frozen phonon at the K-point 
(solid bands) and is compared to that with no frozen phonon 
(dashed bands). Due to the frozen phonon, there are gaps 
at k = which reflect an electron-phonon coupling between 
k = -(2/3)0/a) and k' = (2/3) 0/a). (b) The same calcula- 
tion with the frozen phonon at the F-point. This phonon has 
a coupling which causes forward scattering, (c) The energy 
difference between the electron dispersions with and without 
the frozen phonon at the F-point. The energy difference in the 
nanostrip (solid curve) is compared to that of a (25,25) arm- 
chair nanotube (dashed curve). At low electron wavevectors 
the energy differences are similar, while at high wavevectors 
the energy difference in the nanostrip is suppressed by its edge 
states. 



been folded back to k — 0. The gap is proportional to the 
electron-phonon coupling responsible for backscattering. 
Although it does not reduce the electron-phonon mean 
free path, it is also worth looking at the TO phonon 
at the r-point. This phonon, which only causes for- 
ward scattering, introduces shifts in the dispersion as 
illustrated in Fig.fJJb). Fig.fJJc) shows the energy differ- 
ence between the conduction bands calculated with and 
without the frozen TO phonon. The comparison to an 
armchair nanotube shows that the edge states are essen- 
tially unaffected by the introduction of a frozen phonon 
in the lattice. This result is expected since a wavefunc- 
tion of an edge state is dramatically different from that 
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and the interband matrix elements 
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where L is the length of the strip, m c is the mass of a 
carbon atom, and (ual, uat,ubl, ubt) T is the polariza- 
tion vector normalized to 2. The overall sign depends on 
whether |fe| is smaller or greater than (2/3)(7r/a) but is 
irrelevant for our scattering calculations. 



III. SCATTERING RATES AND MEAN FREE 
PATH 

The electron-phonon scattering rate for an electron 
with wavevector k can be calculated using the Fermi 
Golden Rule 



Y _ p(ab) , p(cm) _ 
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where the absorption and emission terms are 



-,(ab) 
fikq 



2vr, 



h \(k + q,n»-l\H el _ ph \k,n%)\ 2 
x 6(E k+q - E k - TwJ%)p(E k+g ) dE k+q 



and 



i(em) 
flkq 



2tt, 



■\(k + q,n% + l\H el . ph \k,nt q )\ 2 
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where the unoccupied backscattering density of states is 
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where p. is the chemical potential. The phonon occupa- 
tion is given by the Bose-Einstcin distribution 



1 



(11) 



From Eqs. ([5|)- (fT2|) we find that the scattering rate is in- 
versely proportional to the width of the nanostrip. 

Backscattering rates of an undoped (p = 0) GNS at 
room temperature are shown in Fig.fSJa) as a function 
of initial electron energy. Due to the large wavevector 
separation between the initial and final electron states 
the absorbed phonons must have considerable energies 
to ensure energy conservation. These energies are sim- 
ilar to the threshold energies for LA and TO phonon 
emission. The lowest threshold energy corresponds to 
a thermal energy of about T c = huj^ A /k B « 1800 K. 
Therefore, phonon absorption is in a low-temperature 
regime (T <C T c ) where the phonon occupation is ex- 
ponentially small; exp(— Tiu 1 ^ / k-gT) . Because the 
absorption scattering rates scale with phonon occupation 
[see Eq. ©] the backscattering due to phonon absorption 
is expected to be small, even at temperatures well above 
room temperature. This effect does not occur in car- 
bon nanotubes where low-energy long-wavelength acous- 
tic phonons cause backscattering. Backscattering due to 
phonon emission, on the other hand, is of the same order 
as in carbon nanotubes. 

The electron-phonon mean free path is estimated by 
A = |vi|/r where Vi is the initial electron group veloc- 
ity. The mean free path at room temperature is shown 
in Fig.[^b) for a few different chemical potentials. At 
very small energies, the electrons move slowly which ef- 
fectively shortens the mean free path. As the electron en- 
ergy increases, the electron speed approaches the Fermi 
velocity in graphene, vp. Together with increasing inter- 
band LO emission scattering, the mean free path forms 
a maximum. This interband scattering can be reduced 
by doping the material to reduce the number of available 
states in the opposite band. For p — 0.15 eV, the in- 
terband scattering becomes small compared to intraband 
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FIG. 5: (a) Scattering rates of a N zz = 50 zigzag-edge GNS 
at 300 K as a function of electron energy. The LA and TO 
phonons cause intraband and the LO phonon interband scat- 
tering. Intraband phonon absorption occurs at all electron 
energies while intraband phonon emission requires energies 
larger than the threshold energy for the particular branch. 
The threshold energies are about 0.16 eV and 0.18 eV for the 
LA and TO branches, respectively, (b) The total mean free 
path at room temperature due to electron-phonon scattering. 
The mean free path increases with the chemical potential since 
a higher chemical potential reduces the number of unoccupied 
electron states. 



scattering for all energies. The maximum mean free path 
can now be estimated by phonon absorption alone; 



A( ab )(AUT) = A ab) 



N zz 



(12) 



where A, 



(ab) 



70/xm, Ni° z } = 50, and T = 300 K. Above 



the threshold for phonon emission the mean free path 
A( em > is of the order 10 nm which is comparable with 
that in carbon nanotubes^ 2 -. 



IV. CONCLUSIONS 

Using the graphene sheet model [Eq. (JTJ], we have 
found within our model that electron-phonon scattering 
rates in zigzag-edge graphene nanostips is extraordinar- 
ily small for kinetic energies smaller than the phonon 
emission thresholds. The reason for the small scatter- 
ing rates is that long-wavelength acoustic phonons can- 
not cause backscattering of extended states in the lowest 
band. Coupling to the edge states is expected to be small 
since overlaps between extended and edge wavefunctions 
are small and the phonons would have to have a highly 
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exotic symmetry. We have also shown that the electron- 
phonon mean free path is proportional to the width of the 
nanostrip as long as excited bands cannot energetically 
be reached. Wider nanostrips, on the other hand, have 
the disadvantage that they reduce the single-channel win- 
dow AE ps (25/N zz ) eV, and for sufficently wide nano- 
strips (AE ~ k-gT), electrons scattering to the next 
excited band will be appreciable, thereby reducing the 
mean free path. The methodology presented in this paper 
can also be used to calculate electron-phonon scattering 
in armchair-edge nanostrips. The results are expected to 
be similar to zigzag nanotubes, although they might be 
irrelevant as armchair-edge nanostrips are highly sensi- 
tive to edge disorder—. 

Because the electron-phonon mean free path within the 
model is found to be extraordinary long in zigzag-edge 



nanostrips, the actual electron mean free path is likely 
going to be limited either by disorder in the material, de- 
spite a recent prediction that zigzag-edge nanostrips are 
relatively resistant to short-range, long-range, and edge 
disorders^, or by spin polarization effects^ at the edges 
which might open up an additional channel for backscat- 
tering. 
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